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Abstract 

We discuss the extension of the field theoretical approach, already used in the lattice 
determination of the topological susceptibility, to the computation of further terms in the 
expansion of the ground state energy F(9) around 9 = in SU(N) Yang-Mills theories. 
In particular we determine the fourth order term in the expansion for SU (3) pure gauge 
theory and compare our results with previous cooling determinations. In the last part of 
the paper we make some considerations about the nature of the ultraviolet fluctuations 
responsible for the renormalization of the lattice topological charge correlation functions; 
in particular we propose and test an ansatz which leads to improved estimates of the 
fourth and higher order terms in the expansion of F{9). 



(PACS Numbers: 11.15.Ha, ll.10.Gh, 12.38.Gc) 



1. Introduction 



The dependence of SU(N) Yang-Mills theories on the 9 angle is the subject of ongoing 
theoretical debate. The dependence can be expressed in terms of the free energy density 
F{9) which, in the Euclidean theory, is defined as follows: 

exp[-VF(9)] = Z(9) = j [dA\ e S** c (*) e ldQ , (1.1) 

where V is the four dimensional volume, C(x) = \F£ u (x)F®(x) is the usual Yang-Mills 
langrangian and Q = J d 4 x q{x) is the topological charge, with the topological charge 
density q(x) defined as 
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q(x) = ^^u P aF a ^{x)F a pa {x) = d,K,(x) , (1.2) 
where K^x) is the Chern current 

K, = jf^W^ (dpK ~ \9f abc A h p A<)j . (1.3) 
The coefficients of the Taylor expansion of F(9) around 9 = 0, 

fc=0 

are related to the connected expectation values of the topological charge distribution, 

F (k \0) = ^m\e=0 = -i ki ^, (1-5) 

where (-) c is a short notation meaning the connected expectation value taken at 9 = 0. 

F{9) is a non-trivial function, indeed the topological susceptibility, x — {Q 2 )c/V, is 
expected to be different from zero, to the leading order in 1/N, to solve the so-called 
U(l) problem PUI2]- I n Refe. [Sill] it has been argued that, in the large N limit, F(9) is 
a multibranched function of 9, and in particular 

F(9) = F(0) + | min fe {9 + 2nk) 2 + O (1/N) . (1.6) 

Therefore, for sufficiently small values of 9 (9 < n), F(9) is expected to have an almost 
quadratic dependence on 9, with 0(9 4 ) corrections suppressed by powers of 1/N. 
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Numerical Monte Carlo simulations on the lattice are a natural tool to obtain infor- 
mation from first principles about the dependence of F(9) around 9 = 0. While direct 
numerical simulations of the theory at 9 ^ are not feasible because of the complex 
nature of the action, the Taylor expansion of F{9) around 9 = can be computed, in 
principle up to any given order, by measuring the connected expectation values of the 
topological charge over the ensemble of configurations at 9 = 0, as explicited in Eqs. (jl.4|) 
and (jl.5|) . The topological susceptibility has been already extensively studied on the lat- 
tice (see Refs. jSJ |S] for recent reviews) and further terms in the expansion have been 
recently measured [Zj. 

The lattice study of quantities related to topology requires care. The problem is usually 
related to the fact that the topology of gauge configurations on a discretized space-time 
is, strictly speaking, always trivial, and that the usual lattice definition, given in terms of 
gauge fields as a naive discretization of the continuum topological charge, does not have 
the continuum integer valued spectrum: as a good alternative the fermionic definition, 
which is directly related to the index theorem, or the definition based on the so-called 
cooling method, are used. 

In fact, the topological charge operator and its correlation functions can be defined 
on the lattice with the same rigour as for any other operator of the theory: as for any 
other physical quantity, one has to pay attention when removing the ultraviolet (UV) 
regulator, i.e. when going to the continuum limit, since the appropriate renormalizations 
have to be performed. In the cooling method the UV lattice fluctuations, which are 
responsible for the renormalizations, are removed by a process of local minimization of 
the action. However it is also possible to compute the renormalizations and perform the 
appropriate subtractions. This program, usually known as the "field theoretical method" , 
has been already widely discussed and developed, together with a method for the nu- 
merical non-perturbative determination of the renormalization constants, usually known 
as the "heating method" , in the context of the lattice determination of the topological 
susceptibility |H1 El HOI HH H21 H31 H31 HH1 

The aim of the present paper is that of discussing the extension of the field theoretical 
method (and of the heating method used to compute the renormalizations) to the case 
of higher order correlation functions of the topological charge, in order to study the 9 
dependence of the theory. In Section 2, after a review of the field theoretical method, as 
used for the computation of the topological susceptibility, we will discuss its application 
to the case of higher order correlations and develop a suitable extension of the heating 
method. In Section 3 the case of SU(3) pure gauge theory will be used as a testground for 
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the method developed in Section 2, and we will determine the fourth order contribution 
to F{6) and compare our results with those obtained by the cooling technique [7J. In 
Section 4 we will state and test an ansatz about the nature of the UV lattice fluctuations 
responsible for the renormalizations, which will allow us to simplify the computation of the 
connected correlation functions and to obtain more precise determinations of the fourth 
order contribution to F{6). Finally, in Section 5, we will give our conclusions. 



2. Topological charge correlation functions on the lattice 

In this Section we will discuss how the various moments of the lattice topological charge 
distribution renormalize with respect to the continuum ones, and how the corresponding 
renormalizations can be computed numerically In order to make the discussion clearer, 
we will first review the case of the second moment, i.e. the topological susceptibility. 



2.1. Renormalization of the topological susceptibility 

On the lattice it is possible to define a discretized gauge invariant topological charge 
density operator qh(x), and a related topological charge Ql = ^2 x Ql{x) (with the sum 
extended over all lattice points), with the only requirement that, in the formal (naive) 
continuum limit, 

q L (x) a ~Vq(x)+0(a 6 ) , (2.1) 
where a is the lattice spacing. A possible definition is 

-1 ±4 

= 2%2 W^Tr^^n^^)) , (2.2) 

livpcr=±l 

where Yl^ v (x) is the usual plaquette operator in the [iv plane, e^ upa is the standard Levi- 
Civita tensor for positive directions and is otherwise defined by the rule e pupa = — e(_ M )„ p(T . 

A proper renormalization must be performed when going towards the continuum limit, 
like for any other regularized operator. In spite of the formal limit in Eq. (|2.1j) . the 
discretized topological charge density renormalizes multiplicatively jH|: 

q L (x) = Z{(5)a\(3)q{x) + 0(a 6 ) , (2.3) 
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with a multiplicative renormalization constant Z(/3) which is a finite function of the bare 
coupling (3 = 2N/g%, approaching 1 as (3 — > oo. 

When defining the topological susceptibility, further renormalizations can appear. In- 
deed, already the continuum definition, 

X ee M = J d 4 x (q(x)q(O)) , (2.4) 

involves the product of two operators q(x) at the same point: this contact term is divergent 
and not well defined, so that an appropriate prescription must be assigned. It can be 
shown [13 HHj that the correct prescription, corresponding to the quantity which appears 
in the Taylor expansion of F(6), is the one in which the derivative appearing in the 
definition of q(x), Eq. (jl.2|) . is taken out of the vacuum expectation value: 

X = - J d'xd'y dldl{K,(x)K u (y)) . (2.5) 

The lattice definition of the topological susceptibility 

XL = J2^(x)q L (0)) (2.6) 

X 

is in general not guaranteed to meet the correct continuum prescription for the contact 
term, and this leads to the appearance of additive renormalizations: 

XL = Z(p?a*{p)x + M(p), (2.7) 

where M(f3) describes generically the mixing with all local scalar operators appearing in 
the operator product expansion (OPE) of qi^qiiQ) as x ~ in Eq. ([2.6)1 . including in 
particular the action density and the identity operator. 

The idea behind the numerical technique, known as the heating method, used to com- 
pute the two renormalizations Z(/3) and M(/3), is that the UV fluctuations in qL,(x), which 
are responsible for renormalizations, are decoupled from the background topological signal 
so that, starting from a semiclassical configuration of fixed and well known topological 
content, it is possible, by applying a few thermalization steps (i.e. Monte Carlo updat- 
ing steps at the corresponding value of (3), to thermalize the UV fluctuations without 
altering the background topological content. This is certainly true for high enough f3, 
i.e. approaching the continuum limit, and in practice it turns out to be true in a range 
of the f3 values usually chosen in Monte Carlo simulations of Yang-Mills theory, being 
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also favoured by the fact that topological modes have very large autocorrelation times, as 
compared to any other non-top ological mode. 

It is thus possible to create samples of configurations which have a fixed topological 
content, Q, and the UV fluctuations thermalized: various measurements of topological 
quantities on these samples can give information about the renormalizations. For example 
the expectation value of Ql gives 

(Ql) = Z{/3) Q (2.8) 

from which the value of Z((3) can be inferred, while the expectation value of (Q\) gives 

{Ql) = Z((3) Q 2 + V M((3) , (2.9) 

where by V we intend, from now on, the four dimensional volume measured in adimen- 
sional lattice units. 

To check that UV fluctuations have been thermalized, one looks for plateaux in quan- 
tities like (Ql) or (Q L ) as a function of the heating steps performed: only configurations 
obtained after the plateau has been reached are included in the sample. Special care 
has to be paid to verify that during the heating procedure the background topological 
charge is left unchanged: this is usually done by performing a few cooling steps on a copy 
of the heated configuration and configurations where the background topological content 
has changed are discarded from the sample (THj. Further details about the procedure and 
about the estimate of the systematic errors involved can be found in Ref. [Tlj . 

A sample with Q ~ 1 can be used to measure Z and a sample with Q = (usually 
thermalized around the zero field configuration) can be used to determine M. Cross- 
checks can then be performed, using samples obtained starting from various semiclassical 
configurations with the same or different values of Q, to test the validity of the method. 

Once the renormalizations have been computed and the expectation value \l over the 
equilibrium ensemble has been measured, the physical topological susceptibility \ can be 
extracted, using Eq. (|2.7jl . as 

x^m (210) 



If large renormalization effects are present, i.e. if Z <C 1 and if M brings a good 
fraction of the whole signal in xl, the final determination of \i obtained via Eq. (j2.10|) . 
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can be affected by large error bars. However one can exploit the fact that Z and M both 
depend on the lattice discretization qL,{x) and that infinitely many operators qi,{x) can 
be defined all having the same naive continuum limit, to choose improved operators for 
which the renormalization effects are reduced, thus leading to improved estimates of \. 
This is the idea behind the definition of smeared operators JH] 

$(?) = ^S £ WTr (ng(x)ng(x)) , (2.11) 

/i!/p<T=±l 

where 11^ (x) is the plaquette operator constructed with z-times smeared links uji \x), 
which are defined as 

U<?\x) = U,{x), 

±4 

uf(x) = ii-c)u^ 1 \x) + ^J2 u t 1 \^ut 1) ^ + ^ut 1) (x + p:)\ 

a=±l 

Uf{x) = ^ 1/2 , (2.12) 

(^(xyu^(x)) 

where c is a free parameter which can be tuned to optimize the improvement. These 
operators have been successfully used, up to the second smearing level, to determine x 
at zero and finite temperature both in SU{2) [15j (with c = 0.85) and SU(3) [Ej (with 
c = 0.9) pure gauge theory. We will make reference to them later in this paper, when 
computing the higher order correlation functions of the topological charge. 



2.2. Renormalization of higher order correlation functions 

In order to compute the higher order connected moments of the topological charge 
distribution, {Q n ) c , needed for the Taylor expansion of F(9), it is necessary to first com- 
pute the disconnected correlation functions (Q n )- We will consider in particular the case 
n = 4, for which (Q 4 ) c = (Q 4 ) - 3(Q 2 ) 2 . 

Like in the case of the topological susceptibility, also the definition of (Q 4 ) needs a 
special prescription for the contact terms, in order to correspond to the quantity which 
enters the Taylor expansion of F(8). Starting from the lattice definition, Qi, one can 
define the expectation value 

(Ql)= / d 4 X! . . ■d 4 x A {q L (x 1 )q L (x 2 )qL{x^)q L (x4)) (2.13) 
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and it is clear that, apart from an obvious multiplicative renormalization constant Z 4 (/3), 
(Ql) will be linked to (Q 4 ) by additive renormalizations related to the contact terms 
arising when two or more charge densities in the expectation value integrated in Eq. ()2.13j) 
come to the same point. Our aim is to eventually compute the renormalization constants 
numerically, and in this sense we are not interested in the exact form of the mixing terms, 
but just in understanding if they can be related in a simple way to the correlation functions 
of different order, (Q m ), so that they can be computed using the heating method, as it 
will be clarified below. 

We will assume a simple and quite natural form for the renormalization rule: 

(Qt) = Z(^(Q 4 ) + M^((3)(Q 2 ) + M Afi ((3) , (2.14) 

where M 4i2 (/3) and M 4)0 (/3) are two constants which are independent of the topological 
sector. We will discuss in the following whether Eq. (|2.14J) makes sense from a theoretical 
point of view, while its validity will be fully checked numerically in Section 3. 

The presence of the term Z(f3) 4 (Q 4 ) in Eq. (j2.14j) comes from the multiplicative renor- 
malization of the topological charge density. The question is then whether the additive 
renormalizations coming from contact terms can be put in the form M 4:2 ((3){Q 2 )+M 4fi ((3). 
Contact terms arise when two or more charge densities in the expectation value in 
Eq. (j2.13|) come to the same point, but the discussion of the two charge case will suf- 
fice: the instance with three (or more) charge densities coming to the same point can be 
considered as a special case of the two charge case, so that the related mixing terms are 
already considered when discussing the two charge contact term. 

Let us consider for instance £3 ~ x 4 in Eq. ()2.13)) . The operators of lower dimension 
appearing in the OPE of qb^x^qL^Xi) are the identity operator and the action density 
s(x) oc F^ u (x)F^ u (x). The insertion of the identity operator in the expectation value in 
Eq. (j2.13|) leads to a contribution proportional to 

J d 4 Xid 4 x 2 (q L (x 1 )q L (x 2 ) Id) = (Q 2 L ) = Z 2 {(3)(Q 2 ) + V M{(3) , (2.15) 

which is consistent with the assumption in Eq. ()2.14|) . The insertion of the action density 
operator leads to a contribution proportional to 

J d^d'x^xMxMxMx)) = (Q 2 L S) = (S)(Ql) - ±{Ql) , (2.16) 

where S is the total action of the theory and we have used the relation, which is valid for 
any operator O, (d/d/3)(0) = (S)(0) — (SO), where the derivative is taken as the four 



8 



dimensional volume in lattice units, V, is kept fixed. The first term on the right hand side 
of Eq. (j2~TKjl . (S)(Ql), is proportional to Z 2 ((3)(Q 2 ) + V M(f3) and therefore consistent 
with Eq. ()2.14j) . Using the fact that the topological susceptibility \ is a renormalization 
group invariant quantity, i.e. independent of j3, and that \ — (Q 2 )/iVa 4 ), the second 
term can be written as 

^(Ql) = ^(Z\P)xa\P)V + M(P)V)= X V^(Z\P)a\P)) + V^^ = 

= (Q 2 )a~\P)-^(Z^ym + , (2.17) 

which is again consistent with Eq. (|2.14j) . 

This does not complete the discussion, as operators of higher dimension in the OPE 
of q , L(x3)gz,(x4) could bring corrections in which cannot be expressed in the same simple 
form as in Eq. (j2.14j) . However we will assume those corrections to be negligible, and this 
assumption will be well supported by the numerical data shown in Section 3. 

In order to extract the value of (Q 4 ) from Eq. (|2.14|) . we need to measure the lattice 
expectation value (Qf), to already know (Q 2 ), and to determine the new renormalization 
constants M^ifi) and M 4i o(/?)- The computation of the renormalization constants can 
be performed using a simple extension of the heating method. Indeed, measuring the 
expectation value (Q\) on the ensemble thermalized around a semi-classical configuration 
of charge Q, one obtains 

(Q{) = Z\(3) Q 4 + M 4j2 (/?) Q 2 + M 4i0 (/?) . (2.18) 

It is clear that if we repeat the measurement in 2 sectors with different values of Q 
(for instance Q = 0,1), we obtain 2 different constraints involving the renormalization 
constants which, assuming that Z((3) is already known, allow the determination of M 4i2 (/3) 
and M 4 (/3). If we perform the measurement in more than 2 sectors we have more 
constraints than constants to be determined, and this offers the possibility to perform a 
non-trivial test of Eq. (J2.14)) and of the heating method. 

We close this Section by considering the general case of the n-th order correlation 
function. A natural extension of Eq. ()2.14|) is the following: 

n/2 

(Ql) = Z((3) n (Q n ) + J2 M n,n-2k({3)(Q n - 2h ) , (2.19) 

h=l 

and its validity can be discussed along the same lines as for n — 4. In this case, the 
measurement of on the ensemble thermalized around a semi-classical configuration 
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of charge Q gives 

(Ql) = Z n Q n + M n ^ 2 Q n ~ 2 + . . . + M„ i0 , (2.20) 

so that it is necessary to measure i n a ^ l eas t n/2 different topological sectors to 

determine the n/2 renormalization constants M njn _ 2 (/3), M nin _ 4 (/?), . . . , M nj0 (/3). 



3. Determination of (Q 4 ) c in SU(S) pure gauge theory 

In this Section we present numerical results obtained for SU (3) pure gauge theory. We 
will illustrate a detailed study of the renormalization constants involved in the determina- 
tion of (Q 2 ) and (Q 4 ), using the heating method: we will employ samples of configurations 
thermalized in three different topological sectors, Q = 0, 1,2, and this will enable us to 
perform a non-trivial consistency test of Eq. (|2.14j) and of the heating method itself. We 
will then combine the values of the renormalization constants with the results obtained 
at equilibrium to compute (Q A ) C and thus obtain information about the quartic term in 
the expansion of F{6). In particular we will compute the quantity b 2 = — i2(q?) > wrn °h 
measures the relative weight of quartic to quadratic terms and has been determined in 
Ref. [Zj via the cooling method*. 

All the results reported in this Section refer to simulations performed on a lattice 
of size 16 4 , with /3 — 6.1 and the standard Wilson action. Two different discretized 
topological charge density operators have been used, corresponding to the 1-smeared and 
2-smeared operators defined in Section 2. 

3.1. Determination of the renormalization constants 

We will make use of the method described in Section 2 to determine the renormaliza- 
tion constants which enters the computation of (Q 2 ) and (Q 4 )- 

We have collected five different samples of configurations, one thermalized in the Q = 
sector (around the zero field configuration), two in the Q = 1 sector (thermalized around 
two different semiclassical configurations of topological charge one) and two in the Q = 2 
sector (thermalized around two different semiclassical configurations of topological charge 
*We use the same notation used in Ref. [7], in order to make the comparison easier. 
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two). The semiclassical configurations have been obtained by extracting thermalized 
configurations with non-trivial topology from the equilibrium ensemble at /3 — 6.1 and 
then minimizing their action by a usual cooling technique. All the five samples have 
been obtained by performing about 3000 heating trajectories around the semiclassical 
configurations, each trajectory consisting of 90 heating steps; 6 straight cooling steps 
have been applied on heated configurations to check that their background topological 
content did not change. 

We have then measured the expectation values (Qj), and also (Ql)/Q where 

Q / 0, over the five samples, with the aim of applying Eqs. ()2.8j) . ()2.9j) . ()2.18|) and 
determine the renormalization constants*. We have reported the results in Table I for 
the 1-smeared operator and in Table II for the 2-smeared operator: expectation values 
obtained on samples with the same Q turned out to be equal within errors, as they should, 
and we have reported in the tables only their weighted averages. It is still possible to 
appreciate in Table I and II the agreement for the values of Z determined in the Q = 1 
and Q = 2 sectors. We have also reported results for (Q 6 L ), which will be used in Section 
4. 

Let us rewrite Eq. (jZ2J) in the form (Q\) = Z{(3f{Q 2 ) + M 3fl {/3), with M 2fi {(3) = 
VM{(3). The information contained in (Q\) for each different topological sector Q can 
be used to determine M 2i o, using Eq. ()2.9|) . We have 2 constants to be determined, Z and 
M 2j o, and 3 equations (Q = 0, 1, 2), plus two direct determinations of Z from (Ql)/Q, so 
that there are 5 constraints to be satisfied and only 2 variables to be determined. The fact 
that this can be done consistently is a non-trivial test of the method, and more precisely of 
the fact that the UV fluctuations which are responsible for the renormalization constants 
are decoupled from the background topological content, an assumption that is at the very 
basis of the heating method, as it has been explained in Section 2. In practice we have 
used the 5 values measured for (Qj) and (Ql)/Q to perform a best fit to Eqs. ()2.8j) and 
(j2.9|) . obtaining best fit values which are reported in Table III and good x 2 /d.o.f. values 
(~ 0.1 for the 1-smeared operator and ~ 0.2 for the 2-smeared operator*). The values 

^The initial charge of the semiclassical configuration used in the heating procedure, Q, is never strictly 
an integer (apart from the case Q = 0). The deviation from the corresponding integer value is always 
around 3% in our case (e.g., when Q ~ 1, we actually have Q ~ 0.97). When applying Eqs. Ij2.8|l . (|2.9() 
and l|2.18(l . one has to be careful and make use of the real value of Q instead of the closest integer. For 
a more accurate discussion on this point see for instance Ref. |2U| . 

*The low values obtained for x 2 /d.o.f. can be related to the fact that the measurements of {Q\) and 
{Ql)/Q at corresponding values of Q are partially correlated, since they are measured on the same sample 
of configurations. 
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obtained for Z and M 2j0 are compatible with those reported in Ref. [T%] . 

The same procedure has been repeated, using the (Qf) measurements and Eq. (j2.18|) . 
to obtain the best fit values for M 4j2 and M 4j0 reported in Table III. Also in this case we 
have obtained good values for x 2 /d.o.f. (~ 0.1 for both the 1-smeared and the 2-smeared 
operator). This represents a strong numerical support to the validity of Eq. (|2.14j) . 



3.2. Determination of (Q 2 ), (Q 4 ), and (Q 4 ) c 

Now that we have determined the renormalization constants we can proceed to de- 
termine (Q 2 ) and (Q 4 )- The equilibrium values for (Q 2 L ) and (Q\), which are reported 
in Table IV, have been measured on a sample of 300K configurations separated by five 
updating cycles, each composed of a mixture of 4 over-relaxation + 1 heat-bath updating 
sweeps; the reported errors have been estimated by a standard blocking technique. 

The value of (Q 2 ) can be computed as 



(Ql) ~ M. 



(Q 2 ) = ^L—^M . (3.1) 

The expression for (Q 4 ) follows from Eq. (|2.14|) : 

(Qi) = (Ql) - jMg) - M,o (3 2) 

The results for (Q 2 ) and (Q 4 ) are reported in Table IV, the errors have been computed 
by standard error propagation. It is interesting to notice that the values obtained for 
the 1-smeared operator and for the 2-smeared operator are in good agreement, as they 
should, again confirming the robustness of the method. We can finally determine {Q 4 ) c = 
(Q 4 )-3(Q 2 ) 2 , obtaining (Q 4 ) c = 0.32±1.80 for the 1-smeared and (Q 4 ) c = 0.66 ± 0.90 for 
the 2-smeared operator, leading to 6 2 = —0.012(62) and 6 2 — —0.024(32) for the 1-smeared 
and 2-smeared operator respectively, in agreement with the determination reported in 
Ref. 0. 



4. A closer look into the renormalization effects 

The renormalization constants Z, M n>m (m < n) which, for a given lattice discretiza- 
tion Q L , appear in Eq. (|2.19jl . are in principle independent of each other, or at least no 
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simple relation exists among them, unless some further hypothesis can be done about 
the nature of the UV fluctuations which are responsible for the renormalizations. In this 
Section we will propose and test an ansatz which will greatly simplify the structure of 
the renormalization constants and will lead to a renormalization formula which directly 
involves the connected correlation functions, thus allowing a more precise determination 
of b 2 . 

An hypothesis about the nature of the UV fluctuations has been done in Refs. [01 HI] , 
where it was assumed that the discretized topological charge density can be expressed as 

q L (x) ~ [Z + ((x)]q(x) + ri(x) , (4.1) 

where q(x) is a background topological charge density which is determined by physical 
fluctuations on the scale of the correlation length £, whereas ((x) and rj(x) are random 
variables with zero averages which are determined by the short range UV fluctuations and, 
at least in the continuum limit, are expected to be decoupled from q(x), i.e. {((x)q(x)) = 
(r)(x)q(x)) = 0. Summing Eq. ()4.1j) over all lattice points, the following relation follows 
for the lattice topological charge Ql'. 

Q L = ZQ + J2C(x)q(x)+r ] , (4.2) 

X 

where r\ = ^2 x T](x). We now make the further assumption that the term ^2 x C(x)q(x) in 
Eq. ()4.2|) can be neglected, configuration by configuration. This is not unreasonable, in 
view of the fact that q(x) and ((x) are decoupled from each other. We will thus assume 
that 

Q L = ZQ + V , (4.3) 

where rj is a random noise with zero average which is stochastically independent of Q. 

This assumption has relevant consequences for the structure of the renormalization 
constants. Indeed, using the hypothesis that Q and 77 are stochastically independent 
variables and that they are both evenly distributed around zero, it is easy to verify that 
the general renormalization formula in Eq. ()2.19j) becomes 

n/2 , v 

(Ql) = z2[2h) zn ~ 2h{Qn ~ 2h){r,2h) ' (44) 

h=0 ^ ' 

so that the renormalization relation for i s described only in terms of Z and of the 
correlation functions of the noise 77. In particular we have M„ im = (^J Z m (rj n ~ m ) , a 
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relation that should be verified on numerical data if our ansatz in Eq. (J4.3j) is correct. 
From the data in Table III it can be checked that indeed M 42 = 6Z 2 (i] 2 ) = 6Z 2 M 2 fl, 
but we will now proceed further and check the validity of Eq. (|4.4|) up to n = 6. The 
correlation functions of rj can be determined by the heating method using the analogous 
of Eq. (|2.20|) . which up to n = 6 reads: 

(Ql) = Z 2 Q 2 + ( V 2 ) 

(Qi) = Z 4 Q 4 + QZ 2 Q 2 ( V 2 ) + ( V 4 ) 

(Ql) = Z 6 Q 6 + 15Z 4 Q 4 (r ] 2 ) + 15Z 2 Q 2 (r ] 4 ) + (r ] 6 ) . (4.5) 

Using the values for (Ql), (Qj) and (Q L ) obtained in the sectors with Q — 0,1,2 and 
reported in Tables I and II, we have performed a best fit to Eqs. (|4.5jl . obtaining the 
best fit values reported in Table V with x 2 /d.o.f. ~ 0.34 for the 1-smeared operator and 
% 2 /d.o.f. ~ 0.23 for 2-smeared operator. The fact that the values for (Ql), (Qj) and 
(Qi) obtained in the various sectors can be fitted by the simple relations in Eq. (|4.5|) is 
a confirmation of the validity of the ansatz in Eq. (J4.3|) . 

Assuming that Eq. ()4.3|) is valid, it is possible to write a renormalization relation 
which involves directly the connected correlation functions. Indeed, it is a general rule 
that the connected correlation functions of a stochastic variable (Ql in our case), which 
is the sum of two variables which are stochastically independent of each other (ZQ and 
rj in our case), are the sum of the corresponding connected correlation functions, i.e. 

(Q n L ) c = Z n (Q n ) c +( v n ) c . (4.6) 

Therefore in order to compute (Q n ) c we need to know, apart from Z, only one renor- 
malization constant, (f] n ) c , which can be easily measured by computing (Qi) c on the 
sample of configurations in the Q = sector. (Q n ) c is then given by 

(Q n ) c = iQl)c ~ n {Vr% , (4.7) 

where (Q£) c is measured on the ensemble of configurations at equilibrium. 

We have computed (Q\) c = (Qf) ~ 3(<5|) 2 on our equilibrium configurations at (3 = 
6.1, obtaining (Q L ) C = 0.026(7) for the 1-smeared operator and (Q L ) C = 0.057(13) for 
the 2-smeared operator. We have then computed (Q\) c on our sample of configurations 
thermalized in the Q = sector at (3 = 6.1 , obtaining (?7 4 ) c = (t] 4 ) —3(i] 2 ) 2 = 0.006(4) for 
the 1-smeared operator and (f] 4 ) c = 0.001(2) for the 2-smeared operator. In both cases 
(equilibrium and Q = 0) errors have been estimated by standard jackknife techniques. 
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By using Eq. (J4.7j) and the values for Z and (Q 2 ) previously obtained, we have obtained 
(Q 4 ) c = 0.68(24), b 2 = -0.024(9) for the 1-smeared operator and (Q 4 ) c = 0.66(15), 
b 2 = —0.024(6) for the 2-smeared operator. 

By making use of the ansatz in Eq. (J4.3)) we have thus made determinations which 
are much more precise than those obtained in Section 3. The reason is that Eq. ()4.6|) 
allows to relate (Q n ) c directly to the connected correlation functions of the discretized 
lattice topological charge, with only two renormalization constants involved: this greatly 
simplifies computations and error propagation, thus leading to improved estimates. We 
notice that most of the error in the final determination of (Q A ) C and b 2 comes from 
the determination of (Q\) c at equilibrium, which is also the most expensive part of the 
computation in terms of CPU time. The renormalization procedure is thus completely 
under control and numerically non expensive. 

We have also made a determination of b 2 at (3 = 6.0, again on a 16 4 lattice. On a sample 
of about 300K configurations obtained at equilibrium and using the same algorithm as 
for (3 = 6.1 we have obtained, for the 2-smeared operator, {Q 2 L ) = 1.377(7), (Qi) c = 

0. 052(23). On a sample of configurations thermalized in the Q = topological sector 
by performing about 3000 heating trajectories, each composed of 90 heating steps, we 
have obtained, for the 2-smeared operator, (rf) = 0.308(10) and (^ 4 ) c = 0.002(3). From 
these data, using the value Z(f3 = 6.0) = 0.51(2) reported for the 2-smeared operator 
in Ref. [Hj, we obtain b 2 = —0.015(8), which is consistent with the value obtained at 
P = 6.1. 

Let us close this Section with some further speculations about the nature of the UV 
fluctuations. The value obtained for (^ 4 ) c is very small and compatible with zero for 
both the 1-smeared and the 2-smeared operator. We have also measured (?? 6 ) c on the 
sample of configurations at Q = obtaining (i] 6 ) c = 0.001(8) for the 1-smeared and 
(v 6 )c = 0.0005(14) for 2-smeared operator (/3 = 6.1). This seems to indicate that rj 
behaves as a Gaussian noise: it can also be verified on the values reported in Table V that 
(rj 4 ) is always compatible with 3(i] 2 ) 2 and that (rj 6 ) is always compatible with 15(^ 2 ) 3 , as 
it should be for a Gaussian variable. We have also directly verified on the histograms of 
the distribution of t] in the Q = sector that it does not present any significant deviation 
from a Gaussian distribution. If the hypothesis of Gaussian distribution for rj were true, 
the renormalization relation in Eq. ()4.6j) would become, for n > 2, (Q2)c = Z n (Q n )c 

1. e. the connected correlation functions of Ql would renormalize only multiplicatively for 
n > 2. However we will not consider further analysis of this hypothesis in the present 
paper. 
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5. Conclusions 



In this paper we have discussed how to extend the field theoretical method, already 
used for the lattice determination of the topological susceptibility, in order to compute 
further terms in the expansion of the ground state energy F(6) around 6 = 0. 

After a review of the method as used for the determination of the topological suscepti- 
bility, we have discussed the structure of the renormalizations involved in the general case 
and how they can be computed using the heating method. We have presented numerical 
results regarding SU (3) pure gauge theory, providing both support to the correctness of 
the method and a determination of the fourth order term in the expansion of F(6) around 
= 0. 

In the last part of the paper we have made some speculations about the nature of 
the lattice UV fluctuations which are responsible for the renormalizations and proposed 
an ansatz, which we have verified against the numerical data, that leads to a simpler 
structure for the renormalizations and to an improved estimate of the fourth order term. 

The determinations obtained are in agreement with those obtained in Ref. [Zj via the 
cooling method, and confirm that fourth order corrections to the simple 8 2 behaviour of 
F{9) around 9 = are small already for N = 3. 
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TABLE CAPTIONS 



Tab. I. Expectation values measured in different topological sectors for the 1-smeared op- 
erator. 

Tab. II. Expectation values measured in different topological sectors for the 2-smeared op- 
erator. 

Tab. III. Values of the renormalization constants obtained by using the results reported in 
Tables I and II and performing a best fit to Eqs. flUD, O and (EH5j) . 

Tab. IV. Expectation values measured at eequilibrium and results obtained for the renormal- 
ized quantities. 

Tab. V. Values of the renormalization constants obtained by using the results reported in 
Tables I and II and performing a best fit to Eqs. ()2.8|) and ()4.5|) . 
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Table I 



Q 


Z = (Ql)/Q 


{QD 


{QD 


{QD 







0.311(12) 


0.290(20) 


0.48(5) 


1 


0.416(6) 


0.4785(60) 


0.630(15) 


1.295(45) 


2 


0.413(5) 


0.9626(80) 


1.973(50) 


5.82(18) 



Table II 



Q 


Z={Ql)/Q 


{QD 


{QD 


{QD 







0.208(10) 


0.124(10) 


0.131(18) 


l 


0.544(5) 


0.489(5) 


0.556(12) 


0.93(3) 


2 


0.542(4) 


1.314(8) 


2.77(6) 


7.65(17) 



Table III 





operator 


Z 


M 2>0 


M 4>2 


M 4i0 






1 — smeared 


0.414(4) 


0.315(6) 


0.336(16) 


0.289(16) 






2 — smeared 


0.543(5) 


0.211(5) 


0.377(15) 


0.124(9) 












Table IV 












operator 


{QD 


{QD 


{Q 2 ) 


(Q 4 ) 


1 — smeared 


0.7121(38) 


1.548(18) 


2.312(72) 


16.4 ± 1.8 


2 — smeared 


0.8776(60) 


2.368(36) 


2.262(41) 


16.02(72) 



Table V 



operator 


Z 


{v 2 ) 




{v 6 ) 


1 — smeared 


0.415(4) 


0.315(6) 


0.298(11) 


0.462(37) 


2 — smeared 


0.542(4) 


0.211(5) 


0.129(6) 


0.131(17) 
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